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Rotating ground states of a one-dimensional spin-polarized gas of fermionic atoms 
with attractive p-wave interactions on a mesoscopic ring 
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The major finding of this paper is that a one-dimensional spin-polarized gas comprised of an 
even number of fermionic atoms interacting via attractive p-wave interactions and confined to a 
mesoscopic ring has a degenerate pair of ground states that are oppositely rotating. In any realization 
the gas will thus spontaneously rotate one way or the other in spite of the fact that there is no 
external rotation or bias fields. Our goal is to show that this counter-intuitive finding is a natural 
consequence of the combined effects of quantum statistics, ring topology, and exchange interactions. 
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The rapidly increasing sophistication of experimental 
techniques for probing ultracold gases has caused a shift 
of emphasis in theoretical and experimental work in re- 
cent years, from effective field approaches to more re- 
fined methods capable of dealing with strong correlations. 
Such strong correlations arise, for example, in ultracold 
gases confined in dc Broglie waveguides with transverse 
trapping so tight that the atomic dynamics is essentially 
one-dimensional (ID) [Ij, with confinement-induced res- 
onances [l|, y] allowing Feshbach resonance tuning [3\ of 
the effective ID interactions to very large values. This 
has allowed for the experimental verification [J, |5|, |6| of 
the fermionization of bosonic ultracold vapors in such 
geometries predicted by the Fcrmi-Bosc (FB) mapping 
method [7|. The "fermionic Tonks-Girardeau" (FTG) gas 
[8|,|9|, a spin- aligned Fermi gas with very strong attractive 
ID odd- wave interactions, can be realized by a 3D p-wave 
Feshbach resonance as, e.g^ in ultracold "^"K vapor [lOf. 
It has been pointed out [2, Isl, Isj that the generalized 
FB mapping [2, [sj, [oj, lll| can be exploited in the oppo- 
site direction to map the ideal FTG gas with infinitely 
strong attractive p-wave interactions to the ideal Bose 
gas, leading to "bosonization" of many properties of this 
Fermi system. One of us (MDG) in collaboration with 
A. Minguzzi showed recently il2j that for an even num- 
ber A^ of fcrmions on a mesoscopic ring, the ideal FTG 
ground state is highly degenerate due to fragmentation 
of the mapped ideal Bose gas between two macroscopi- 
cally occupied orbitals with circumferential wave vectors 
±7r/L, where L is the ring circumference. 

The purpose of the present paper is to examine how 
the ground state properties of a FTG gas on a mesoscopic 
ring are changed for an even number N of fermions when 
the strength of the atom-atom attraction is made finite, 
although still very large, in which case there exists an 
FB mapping between the FTG gas and a system of N 
bosons with repulsive delta- function interactions. Based 
on this we have found the remarkable result that whilst 
most of the ground state degeneracy of the ideal FTG 
gas is lifted a twofold degeneracy remains, and these two 



states have nonzero total angular momentum zL^Nh, cor- 
responding to BEG of all N particles of the mapped Bose 
gas into either an orbital of angular momentum ^ or one 
of angular momentum — ^. The lifting of the degener- 
acy will be given a natural interpretation in terms of the 
mapped Bose system using the famous no-fragmentation 
theorem of Nozieres and Saint James [l3|. Thus, in any 
given experimental realization with an even number N 
of fermions the FTG will spontaneously rotate one way 
or the other around the ring with equal probability. This 
apparently violates a theorem of F. Bloch Jdj requiring 
the ground state to be non-rotating, and we shall demon- 
strate why the theorem fails for our case of fermions on 
a mesoscopic ring. 

Ideal FTG gas: To set the stage we consider the ideal 
FTG gas in general and on a ring. The ideal FTG gas 
consists of spin-aligned fermions with infinitely-strong 
zero-range attractions which are a zero-width, infinite 
depth limit of a square well of depth Vq and width 2a;o, 
with the limit taken such that VqXq — > (Trh?)/8meff 
where rueff is the effective mass of the colliding pair 
@j@,H- It causes odd- wave scattering (ID analog of 3D 
p-wave scattering) with ID scattering length aiD = — oo, 
with the result that all energy eigenstates ^f of the FTG 
gas are obtained from corresponding ideal Bose gas states 
^B by the mapping ^s -^ "ip = A^^b, where A is the 
same mapping function used in the original solution of 
the TG gas JTl: A{xi, ■ ■ ■ ,xn) = ni<j<c<Ar sgn(a:j -xe) 
where the sign function sgn(x) is +1 (—1) if a; > (a; < 
0). This introduces sign-changing discontinuities in the 
FTG states necessary to reconcile fermionic antisym- 
metry with a strong interaction in the zero-range limit 
Xq ^ [llj. Inside the square well the solution passes 
smoothly through a zero at Xj — xg — 0, so the discon- 



tinuity is an illusion of the zero-range limit |2|. I8l.l9llll|. 
The exact ground state is 
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where (j)o{x) is the lowest single-particle orbital for the 
given boundary conditions, illustrating the mapping from 
the ideal Bose gas ground state "^bo = Hj=i't'o{xj) 
to the fermionic FTG ground state "ipo- The physical 
meaning of this mapping is clarified by looking at the so- 
lution of the two-fermion problem both inside and outside 
the square well, before passing to the limit xo -^ 0. De- 
noting the relative coordinate by X12 = a;i — 2:2, the exte- 
rior solution (1x12 1 > xo) in the ideal FTG limit (ID scat- 
tering length aio — > —00) represents a zero-energy scat- 
tering resonance, with wave function -1-1 for a;i2 > and - 
1 for a;i2 < 0, and the interior solution (|a;i2 1 < 2:0) fitting 
smoothly onto this is sin(Ka;i2) with k = y^mVo/h? — 
7r/2a;o. The corresponding mapped Bose gas state is 
-1-1 everywhere outside the well (ideal Bose gas ground 
state), but inside the well it is sin(K|a;i2|). Since this 
vanishes with a cusp at a:i2 = 0, physical consistency 
requires the presence of a zero-diameter hard core inter- 
action added to the square well. The mapped Bose gas 
is then not truly ideal, but rather a TG gas with super- 
imposed attractive well, whose nontrivial interior wave 
function becomes invisible in the zero-range limit, simu- 
lating an ideal Bose gas insofar as the energy and exterior 
wave function are concerned. The required impenetrable 
core is physically quite reasonable, since the atoms have 
a strong short-range Pauli exclusion repulsion of their in- 
ner shells, whose diameter is effectively zero and strength 
infinite on length and energy scales appropriate to ultra- 
cold gas experiments. 

Ideal FTG gas on a ring: Here we review the discus- 
sion of Girardeau and Minguzzi for fermions on a ring of 
circumference L [12l | . For an odd number A'^ of fermions 
Eq. ^ is the correct FTG ground state, with 0o = 1/V^ 
the trivial and unique ground state orbital, a plane- wave 
with zero wave vector k and hence zero angular momen- 
tum. However, if N is even we encounter a difhculty. 
For example, for N = 2 and < X2 < L the mapping 
function A = sgn(a;i — X2) is -1 for zi = -I- e but -1-1 
for xi = L — e, where e — > 0+. As a result, ^_fo is 
not periodic but instead antiperiodic, violating the re- 
quirement that the fermion wave function be single val- 
ued. More generally, A{xi,--- ,xn) is antiperiodic for 
all even N, leading to the same difficulty. This problem 
is easily repaired by taking the Bose ground state ^^o 
to be the lowest ari^jperiodic state, a fragmented BEG 
with wN atoms in the orbital e^^^i/^ and (1 — w)N in 
^-i-RXj/L r^[^Yi < w < 1. The fermionic ground state 
"ifpo — A'^BQ is then properly periodic, but is {N + 1)- 
fold degenerate and conveniently labeled by a quantum 
number i, = {w - ^)N = 0,±1,±2,--- ,±f related 
to the eigenvalue Lz of orbital angular momentum by 
Lz = izfi- Since N is even, Lz takes on all integral values 
from ~^Nh to ^Nh, in spite of the half-integral angular 
momenta of the two orbitals, and the fermionic ground 
states have the same angular momentum spectrum since 
the total angular momentum operator commutes with 



the mapping function A. 

The above discussion leading to rotating and degener- 
ate ground states clearly relies on the mesoscopic nature 
of the ring as it relies on the discreteness of the angular 
momentum. In contrast, Bloch's "no-rotation" theorem 
14| is based on a Galilean transformation treating the 
momentum as a continuous variable and therefore holds 
strictly only in the thermodynamic limit; that is why 
our mesoscopic system can violate Bloch's theorem. Fur- 
thermore, the physics underlying why FTG gases with 
an even number of fermions on a ring can have rotating 
ground states is also revealed in the above discussion: 
It is a consequence of the combined effects of quantum 
statistics, which is implicit in the use of the mapping 
function A{xi , • ■ • , xjf) to incorporate the anti-symmetry 
of the fermion wave function, and the ring topology which 
necessitates a non-zero angular momentum so that the 
underlying Bose ground state can exhibit the antiperiodic 
boundary conditions needed to counter the antipcriodic- 
ity of A. The same general argument applies also for the 
case of finite interactions since the problem can still be 
approached using the same mapping function A (see be- 
low), albeit that the mapped Bose problem is no longer 
ideal. However, for the ideal FTG it is still possible to 
consider combinations of the degenerate ground states 
that have zero angular momentum and to address this 
we need to consider the effects of finite interactions. 

FTG with finite interactions: It has been realized in re- 
cent years that the FB mapping method used to exactly 
solve the TG gas [2] and FTG gas [SJ, |9[ is not restricted 
to the TG case of point hard core boson-boson repul- 
sion and the FTG case of infinite zero-range fermion- 
fcrmion attraction. In fact, exactly the same unit anti- 
symmetric mapping function A(xi,--- ,xn) used there 
also provides an exact mapping between bosons with 
delta-function repulsions g^j-,5{xj — xe) of any strength, 
i.e., the Lieb-Liniger (LL) model [l5|, and spin-aligned 
fermions with attractive interactions of a generalized 
FTG form with reciprocal fermionic coupling constant 
g[ij- The ID scattering length oid < is invariant under 
the mapping, which, therefore, maps strongly attractive 
fermions to weakly repulsive bosons and strongly repul- 
sive bosons to weakly attractive fermions, the coupling 
constants being related to the scattering length and each 
other by g^^ = 2h? /m\aiD\, fffo = 2/i^|ai£)|/m, and 
9?d9id = ^ i i i El- The physical meaning of 
this mapping is illustrated by FiglTJ which shows, for the 
case N — 2, how the mapping converts a FTG attractive 
interaction to the LL delta function interaction, which 
generates cusps in the wave function in accord with the 
LL contact conditions [l5|. The potential, both for ^^ 
and ^B, consists of a deep and narrow square well of 
width 2x0 and depth Vq plus a TG impenetrable core at 
a;i2 = 0, as discussed previously for the ideal FTG limit 
o-iD — *■ —00. The point core has no effect on ^p, which 
already vanishes at a;i2 = due to antisymmetry, but it 
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¥b = sgn(Xi2) v|/p 



V, = sgn(Xi2) Vb 



FIG. 1: Two-particle fermionic wave function ^f and mapped 
bosonic wave function ^s as a function of xi for fixed 12 • The 
potential for both ^f and 4's is a deep and narrow square 
well plus a point hard core, but in the zero-range limit its 
effect on "^b outside the well is the same as that of vb ~ 
gio3ixi2). 



is necessary for ^s for consistency of the mapping. The 
generalized FTG well involves going to the zero width 
limit xo ^ and infinite depth limit I/q ^ 00 with 
VqXq held constant as for the FTG limit auj — > — cx), 
but now at a slightly smaller value, such that the pa- 
rameter K in the interior wave function sin(KXi2) scales 
as K 



7, — —. 1 as a;o ^ fi 01 . Outside the well 

this potential has the same effect on ^s as a LL delta 
function (7^(5(a;i2); this is the physical meaning of the 
mapping from the spin-aligned Fermi gas to the LL Bose 
gas. Here we will use this mapping to treat the case 
of spin-aligned fermions with strong but finite attraction 
via mapping to an LL gas with weak repulsions. 

The relative strengths of the FTG attractive p-wave 
interactions and the LL repulsive interactions are quan- 
tified in the dimensionless coupling coefficients 7f — 



2£idI 



^'f and 7b = -^^^f-, respectively, where n = ^ is 
the linear atomic density, and 7f7b = 4. The ideal FTG 
gas is realized in the limit aiD —* —00 so that 7^? ^ 00 
and 7b — > 0. Here we are interested in the limit of fi- 
nite but still very strong p-wave attractive interactions, 
or the limit 7f >> 1 and 7_b << 1. Moreover, this 
is precisely the limit where the ground state properties 
of the LL Hamiltonian can be accurately captured us- 



ing mean field theory in which all N bosons are assumed 
to occupy the same single particle orbital 4> that is de- 
termined self-consistently via the solution of the Gross- 
Pitaevskii (GP) equation [l6|. The A^-particle boson 
wave function is then written in second quantized form as 
I*i3) = (at)^|0)/VM, with a = ^ dx(j)*{x)i){x), i)^x) 
and 'ijjiyX) being boson creation and annihilation field op- 
erators. Lieb et al. [17| have studied the parameter space 
in which a system of N bosons with repulsive interactions 
in a 3D trap with tight transverse binding may be treated 
as an effective ID LL gas, and the GP approach should be 
valid in the limit jb < l/-^^i the same criterion applying 
to our ring trap with tight transverse binding. 

In order to address fragmentation of the ground state 
of the FTG gas in the presence of finite p-wave inter- 
actions, here we adopt a more general Hartee-Fock ap- 
proach to the mapped LL problem in the limit 75 <C 1, 
which is a two-component generalization of the GP the- 
ory, where the two components are the two degenerate 
ground state orbitals (j)±{x) = e±*''^/^/\/L of the GP 
equation with antiperiodic boundary conditions and en- 



j^N{N - I) where g = .gf^,. 



exact second-quantized LL Hamiltonian is 
H= I dxl-^'iIjUx)-^ipix 



The 
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and we evaluate its expectation value in the (A^-t-l) states 
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where a± — J dx(j)'^{x)^p{x). The states in Eq. ^ rep- 
resent fragmented states of the two orbitals 4>±{x) in 
second-quantized form, and are the analogue of the pre- 
viously discussed fragmented states ^so associated with 
the ideal FTG gas. The energy expectation values of the 
LL Hamiltonian ^ taken with respect to the fragmented 
states Q are given by 
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The energy assumes a maximum at u; = i, represent- 
ing a fragmented state with zero total angular momen- 
tum and equal occupations ^ '^^ the orbitals (f>±, and 
the energy minima are at w = 1 and w = 0, the un- 
fragmented states with all N atoms in either 0+ or (/>_, 

energies -Bo(l) = -Bo(O) = ^N + ^NiN - I), and 
total angular momenta ±-^. These two rotating states 
lie lower than the nonrotating state w = ^ by an amount 
-^N"^ , the extra exchange energy arising from fragmen- 
tation of the state w = ^, in accordance with the no- 
fragmentation theorem [13| . which can be regarded as a 
physical mechanism forcing the ground LL state to be 
rotating. Since the total energy and total angular mo- 
mentum operators commute with the mapping function 



A, it follows that the FTG gas has two degenerate ground 
states with opposite angular momenta ±-^, and this is 
the main result of this paper. 

It is instructive to compare these mean field results 
with exact solutions of the LL Bethe ansatz equations 
[15| for N ~ 2 and antiperiodic boundary conditions. 
These states have the form '^{xi,X2) = e^'^^^^e'^'^^^'^ — 
^le ^ik^xi ^ikix2 jjj ^jjg ggctor Q < xi < X2 < L, and are 

equal to ^(0:2, xi) when Q < X2 < xi < L. For 75 <C 1 
the lowest nonrotating state has ki — ~k and k2 = k 
where kL = n + e with < e ^ 1, and we find e — 
^ + 0(7|). Its energy is 2 (f )^ + ^ + 0(jI), which 
agrees to order ^ with the mean field energy £^0(5) of 
Eq. (gl), with LL units ?i = 2m = I and 5 = 2c = ^. 
The twofold degenerate ground states ^q have kiL ^ 
±7r — e and k2L — ±7r + e with < e ^ 1, and we 
find e = V^Ts + 0(75) and energy 2 (^) + ^, which 
again agrees with the mean field result £^0(0) = £'o(l) of 
Eq. ([U to order ^. We conclude that our mean field 
solutions capture both the qualitative and quantitative 
features of the exact antiperiodic LL solutions for 73 <C 
I, and in particular, the exact LL ground state is twofold 
degenerate and rotating, with angular momenta ±^-- 

A potential scheme to measure the angular momen- 
tum of the rotating gas is the non-destructive approach 
discussed in Ref. [18], which employs Raman transitions 
between hyperfine levels of the atoms due to oppositely 
circularly polarized laser fields that propagate along a 
direction in the plane of the atomic ring. The spatial ab- 
sorption image for the circular polarizations can depend 
on the angular momentum state of the gas since the rota- 
tion gives rise to a Doppler shift that modifies the proba- 
bility of the Raman transition. Then if the gas is rotating 
the absorption image will show an asymmetry with re- 
spect to the axis passing through the center of the ring, 
since the atoms on either side of the axis will be mov- 
ing in opposite directions, experience different Doppler 
shifts, and hence produce different absorptions. 

Ring dark solitons: So far we have tacitly assumed that 
the density profile of the ground state of the FTG gas, 
and the underlying LL gas, is spatially homogeneous. Re- 
laxing this assumption raises the possibility that a non- 
rotating ground state may be realized in the form of a 
ring dark soliton. Dark solitons are real valued solutions 
4>Ds{x) of the GP equation [l9|, neglecting any overall 
phase factor, which have an inhomogeneous density pro- 
file and a phase change of tt at the position where the den- 
sity touches zero. Thus dark solitons have the required 
antiperiodicity for the mapped Bose solutions. Details 
of how to calculate the ring dark solitons for finite inter- 
actions are given exhaustively in Ref. 19], and for the 
ideal FTG the dark soliton solution of the GP equation 

with {g ~ 0) degenerates into 4'Dsix) — \/ x sin(7rx/L), 
and is degenerate with the previous rotating plane-wave 



solutions. When the ring dark solitons {g > 0) arise they 
always have energies higher than the rotating plane- wave 
solutions, so that although they are non-rotating they 
represent excited states. An intuitive way to convey this 
is that the healing length associated with the ring dark 
solitons always engenders a larger kinetic energy penalty 
than the wavevectors ±7r/L associated with the rotating 
ground states. 

In summary, we predict that the ground state of an 
FTG gas with an even number of fermions on a ring 
and finite p-wave interactions is doubly degenerate and 
rotating. We have further shown that the rotation is 
a macroscopic manifestation of the strong microscopic 
many-body correlations present in the gas. Observation 
of the spontaneous rotation of an FTG gas on a ring 
would constitute a striking validation of the underlying 
strong correlations. 



* Electronic address 
t Electronic address 



girardeau@optics.arizona.edu| 



'ewan.wright@optics.arizo na.edu] 
[1] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[2] B.E. Granger and D. Blume, Phys. Rev. Lett. 92, 133202 

(2004). 
[3] J.L. Roberts et al, Phys. Rev. Lett. 86, 4211 (2001). 
[4] B. Paredes, et al, Nature 429, 277 (2004); T. Kinoshita, 
T.R. Wenger, and D.S. Weiss, Science 305, 1125 (2004). 
[5] T. Kinoshita, T.R. Wenger, and D.S. Weiss, Phys. Rev. 

Lett. 95, 190406 (2005). 
[6] T. Kinoshita, T.R. Wenger, and D.S. Weiss, Nature 440, 

900 (2006). 
[7] M. Girardeau, J. Math. Phys. 1, 516 (1960); M.D. Gi- 
rardeau, Phys. Rev. 139, B500 (1965), Sees. 2, 3, and 
6. 
[8] M.D. Girardeau and M. Olshanu, cond-mat/0309396 
[9] M.D. Girardeau, Hieu Nguyen, and M. Olshanii, Optics 
Communications 243, 3 (2004). 
[10] C. Ticknor, C.A. Regal, D.S. Jin, and J.L. Bohn, Phys. 

Rev. A 69, 042712 (2004). 
[11] T. Cheon and T. Shigehara, Phys. Lett. A 243, 111 

(1998) and Phys. Rev. Lett. 82, 2536 (1999). 
[12] M.D. Girardeau and A. Minguzzi, Phys. Rev. Lett. 96, 

080404 (2006). 
[13] P. Nozieres and D. Saint James, J. Phys (Paris) 43, 1133 

(1982). 
[14] D. Bohm, Phys. Rev. 75, 502 (1949). 
[15] E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 
[16] See, for example, E. M. Lifshitz and L. P. Pitaevskii, 
Statistical Physics, Part 2 (Pergamon Press, Oxford, 
1989)pp. 85-118. 
[17] E.H. Lieb, R. Seiringer, and J. Yngvason, Phys. Rev. 
Lett. 91, 150401 (2003); see REGION 2, p. 3 of this 
paper. 
[18] K. P. Marzlin, W. Zhang, and E. M. Wright,Phys. Rev. 

Lett. 79, 4728 (1997). 
[19] L. D. Carr, C. W. Clark, and W. P. Reinhardt, Phys. 
Rev. A 62, 063610 (2000). 



